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ABSTRACT 


V  !zn-U*  >>P.S 

A  Poisson  stream  of  arrival  rate  and  service  time  distribution 

-  •  -f  ,  ,  I 

"  y  / 
Aj(x)  has  preempt  resume  priority  over  a  second  stream  of  rate  Xjj  >  ' 

and  distribution  A^j(x).  Abundant  theoretical  results  exist  for  this 
system,  but  severe  numerical  difficulties  have  made  many  descriptive 
distributions  unavailable.  Moreover,  the  distribution  of  total  time 
in  system  of  low  priority  customers  has  not  been  discussed  theoretically. 
It  is  shown  that  the  waiting  time  sequences  of  such  customers  before 
first  entry  into  service  is  a  Lindley  process  modified  by  replacement. 
This  leads  to  the  total  time  distribution  needed.  A  variety  of  des¬ 
criptive  distributions,  transient  and  stationary,  is  obtained  numer¬ 
ically  via  the  Laguerre  transform  method.  -- 


Key  words:  Priority  queue.  Preempt /resume  discipline,  Modified 

Lindley  process,  total  time  in  system,  Numerical  distri¬ 
butions 
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§0.  Introduction 


A  priority  queue  describes  two  M/G/l  service  systems  interacting 
through  a  common  server.  The  first  system,  to  be  called  system  1,  has 
a  Poisson  input  stream  of  customers  of  intensity  A ^ ,  and  independent 
service  times  with  c.d.f.  A^(x).  Similarly,  system  II  has  intensity 
Ajj  and  c.d.f.  Ajj(x).  All  interarrival  times  and  service  times  are 
independent  of  each  other.  A  substantial  literature  dating  back  to  the 
fifties  treats  the  interaction  of  the  two  customer  streams  when  system  I 
customers  have  preemptive  priority  over  system  II  customers.  Abundant 
results  have  been  obtained  giving  first  and  second  moments  of  interest 
and  some  asymptotic  results  (see,  e.g.,  Gaver  [1],  Heathcote  [5], 

Jaiswal  [4],  Keilson  [5,8,9],  Miller  [18]  and  Prabhu  [19]].  The  distri¬ 
bution  of  random  variates  of  interest  required  for  system  design  have  not 
been  available,  however,  because  of  Laplace  transform  difficulties.  Many 
asymptotic  results  such  as  heavy  traffic  approximations  for  waiting  times 
have  been  flawed  by  lack  of  error  bounds  and  disturbing  relaxation  time 
problems. 

An  algorithmic  procedure  has  been  needed  providing  accurate  numerical 
distributions  for  effective  service  times,  busy  periods,  and  ergodic 
waiting  times.  The  Laguerre  transform  method  introduced  by  Keilson  and 
Nunn  [11],  Keilson,  Nunn  and  Sumita  [12]  and  studied  further  by  Sumita  [21] 
provides  a  framework  for  evaluating  multiple  convolutions  and  other  con¬ 
tinuum  operations.  Many  of  the  distributions  required  have  been  obtained 
previously  thereby  [11,12,13,14,20,21].  For  the  total  time  in  system  of 
low  priority  customers  in  the  system  of  interest,  however,  new  probabilistic 


analysis  has  been  needed  before  the  Laguerre  procedure  is  applicable. 
This  analysis  and  related  asymptotic  results  are  the  focus  of  this 
paper.  The  Laguerre  transform  will  serve  only  as  a  tool,  but  one  whose 
power  will  be  made  evident  hopefully  through  the  results. 

A  recent  paper  by  the  authors,  "The  Depletion  Time  for  M/G/l  Sys¬ 
tems  and  a  Related  Limit  Theorem"  [14],  discusses  single  server  M/G/l 
systems  with  many  classes  of  customers  and  complex  order  of  service  dis¬ 
ciplines.  That  paper  provides  a  survival  function  bound  for  the  time  in 
system  of  any  customer  at  ergodicity,  giving  rise  to  a  robust  (but  non¬ 
exponential)  limit  theoretic  bound  for  heavy  traffic.  In  the  present 
paper,  the  significance  of  the  service  time  distributions  of  competing 
classes  and  of  their  traffic  intensities  is  emphasised. 

In  Section  1,  the  system  studied  is  described  and  the  main  results 
are  summarized.  A  modified  Lindley  process  is  derived  in  Section  2, 
which  represents  the  waiting  time  until  first  entry  into  service  of  the 
k-th  low  priority  customer.  In  Section  5,  the  Laplace  transforms  of 
stationary  distributions  of  interest  are  given,  and  related  heavy 
traffic  limit  theorems  are  established.  A  final  section  is  devoted 
to  numerical  examples.  All  ergodic  and  transient  distributions  des¬ 
cribed  in  the  previous  sections  are  evaluated  numerically  using  the 
Laguerre  transform  method. 


-  - 
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§  1 .  The  system  of  interest  and  main  results 

We  consider  the  two  customer  streams  described  in  Section  0  with 
arrival  intensities  X^  and  X ^  and  service  times  and  with  c.d.f. 's 
Aj(x)  and  Ajj(x).  We  assume  that  Tj.  and  have  finite  second  moments 
and  that  the  ergodicity  condition  p<.  =  +  p^  <  1  holds,  where 

PI  =  anc*  Pjj  =  A 1 1 E  TTi  1 1  *  corresponding  transforms  are 

-wTi  ~W^II 

denoted  by  cij(w)  =  E[e  ]  and  a^(w)  =  E[e  ].  A  class  I  customer 

evicts  any  class  II  customer  from  the  service  facility.  When  that  class 

I  customer  and  all  subsequent  class  I  customers  with  overlapping  presence 

complete  service,  the  evicted  class  II  customer  resumes  service.  The 

queue  discipline  of  class  II  is  FIFO.  The  queue  discipline  in  class  I 

is  irrelevant  for  class  II  service  delay  provided  the  busy  period  for 

class  I  customers  is  undisturbed. 

A  tool  we  shall  appeal  to  frequently  is  that  of  effective  service 

time  distribution.  Suppose  a  random  service  time  T  in  the  absence  of 

interruptions  has  c.d.f.  A(x).  Interruptions  A  with  c.d.f.  B(xl  arrive 

in  a  Poisson  stream  of  rate  X.  Let  a(w)  =  E[e  and  6(w)  =  E[e  wA] . 

eff 

Then  it  has  been  shown  [4,  5]  that  the  elapsed  time  T  until  service 

eff 

eff  eff  -wT 

is  completed  has  c.d.f.  A  (x)  with  a  (w)  =  E[e  ]  given  by 


(1.1)  cteff(w)  =  a (w  +  X  -  X0(w))  . 


In  particular,  for  the  preempt-resume  system  described,  the  effective 

eff 

service  time  transform  (w)  of  class  II  customers  is  given  by 


(1.2) 


eff 

an  (w)  =  an(w  *  xi  '  Vbpi^ 
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since  class  I  customers  ignore  all  class  II  customers  and  see  their  own 
M/G/l  system.  Here  o^pjCw)  is  the  server  busy  period  transform  for  class  I 
M/G/l  systems.  This  transform  agpj  (V)  satisfies  the  classical  Takdcs  equa¬ 
tion  [22] 


(1.3) 


oBpi(w)  =  «j(w  +  -  XI°BPI(W))  • 


Of  related  interest  is  the  stationary  waiting  time  distribution  for 
class  I  and  class  II  customers.  That  for  class  I  is  given  by  the  familiar 
Pollaczek-Khintchine  distribution  Fp^jC*)  [lb]  with  transform 


-wlV 


(1-4) 


♦pkl'"’  *  E[e 


PKI 


1 


1  = 


I 


1  -  aI (w) 

P  t  ^  , .  i '  r  ,T'  t  ] 


I1  wEtTj] 


where  WpKI  is  the  stationary  waiting  time  for  class  I.  For  the  class  II 
customers,  the  discussion  of  waiting  time  is  much  more  difficult.  Such 
a  customer  has  a  waiting  time  before  first  entry  into  service,  and  may  be 
evicted  repeatedly  by  streams  of  overlapping  class  I  customers.  For  the 
class  II  customers,  the  total  time  in  system  rather  than  the  waiting  time 
is  needed.  Nevertheless,  the  waiting  time  before  first  entry  into  service 
provides  a  stepping  stone  to  the  time  in  system.  Let  be  the  time  until 
first  entry  into  service  of  the  k-th  class  II  customer.  It  will  be  shown 
in  Section  2  that  the  sequence  has  a  Lindley  process-like  structure  in 
that  one  has 


(1.5) 


k+1 


K.  + 


k  £k.l  if  Bk  *  ?k.l  2  0  • 


T 

XO 


if  hk  *  sk.i  2  0  • 
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K' ' 


u 
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eff  eff 

Here  £,  ,  =  S.  -  A,  ,  where  S.  is  the  effective  service  time  of  the 
k+1  k  k+1  k 

k-th  class  II  customer,  i.i.d.  with  transform  given  in  (1.2),  and  is 

the  interarrival  time  between  the  k-th  and  the  (k+l)-st  class  II  customers. 

The  variate  T^  is  the  first  passage  time  of  the  server  backlog  process 

Bj(t)  of  class  I  M/G/l  system  from  B^ (0+)  =  X  to  zero  [14].  The  variate 

X  is  given,  as  we  will  see  in  Section  2,  by  X  =  B T  (v^—  E)  with  E  being  the 

I  An 

exponential  variate  of  the  unit  mean.  The  process  may  be  called  a 
Lindley  process  [i“]  modified  by  replacement.  An  algorithmic  procedure 
will  be  given  for  evaluating  the  sequence  of  distributions  of  K,  recursively 
via  (1.5)  based  on  the  Laguerre  transform  method,  and  in  turn  to  the  sta¬ 
tionary  distribution  of  .  The  distributions  of  are  of  separate  inter¬ 
est  in  that  they  display  the  approach  to  ergodicity  and  provides  relaxation 
time  information  numerically. 

Alternatively,  it  will  be  shown  that  the  stationary  waiting  time 
before  first  entry  into  service  of  the  class  II  customers  is  given  by 


(1.6)  <J>n(w)  =  $PKS(W  +  '  xiaBPICw)) 


-wW 


where  <J>pKS(w)  =  E[e 


PKS- 


]  and  Kpt,g  is  the  stationary  waiting  time  seen  by 
all  customers  in  the  system  when  there  is  no  priority  and  the  service  dis¬ 
cipline  is  FIFO.  A  comparison  of  (1.6)  with  (1.1)  then  gives  rise  to  the 
following  formal  statement. 


At  &tcuUonasUty  a  cJLa&A  II  cuAtomeA  zxpeAizncu ,  before 
its  fiA&t  entry  into  service,  the  system  P-K  delay  for 

(1.7)  xs  =  AI  +  AII  and  AsCx)  =  fxiAi(x^  +  ^nAn(x)]/As  modi¬ 
fied  by  interruption &  at  Poisson  note.  A  tilth  duration 
T 

]BPr 
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§2.  The  modified  Lindley  process  with  replacement  for 

In  this  section  we  establish  the  recursion  relations  (1.5)  for  the 
waiting  times  K^  of  the  class  II  customers  described  in  Section  1.  The 
total  time  in  system  of  the  k-th  class  II  customer  is  then  obtained  using 
WR  as  a  stepping  stone. 

For  notational  convenience,  we  denote  the  k-th  class  II  customer  by 
C^.  Let  arrive  at  the  system  at  time  Ke  note  that  the  interarrival 

times 


(2.1) 


k+1 


k+1 


-  r. 


are  i.i.d.  and  exponentially  distributed  with  parameter  Xjj.  Let  first 

*  *  * 
receive  service  at  time  Tj,  and  leave  the  system  at  time  so  that 


(2.2) 


Tk  =  Tk  +  Wk  5 


**  *  ceff  u-  ,  ceff 

T,.  =  t.  +  S,  =  T,  +  K,  +  s,. 


eff 

Here  ,  the  effective  service  time  of  C^,  is  i.i.d.  with  transform 

eff  -wSe^  eff 

ajj  (w)  =  E[e  ]  given  by  (1.2).  Ke  note  that  h'k  and  Sj,  are  inde¬ 

pendent  and  that  the  total  time  in  system,  U^,  of  is  given  by 

<2-5>  uk  *  \  \  -  "k  *  skff  • 

eff 

Since  the  distribution  of  is  known,  the  distribution  of  K^  leads  to 

that  of  U.  . 

k 

Suppose  that  is  , ready  present  when  leaves  the  system  at 

**  eff 

.  This  is  equivalent  to  saying  that  +  ^  s  =  K^  +  .  Hence 

K'k+j  is  given  by  (cf.  Figure  2.1a), 


HO- 
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(2'4)  "k+1  "  "k  +  Ck+1  lf  "k  +  ?k+l  ~  0 


where 


(2-5) 


5k+i  =  sk  -  \+l 


We  note  from  (2.2)  and  (2.4)  that  the  first  service  entry  epoch  of 
when  K  +  >0,  is 


(2.6) 


,  =  T,  ,  +  W,  ,  =  T,  +  W.  +  S,  =  T,  ,  W,  +  £.  ,  >  0  , 

c+1  k+1  k+1  k  k  k  k  k  n+1 


since  no  class  I  customers  are  present  at  departures  of  class  II  customers. 

We  now  suppose  that  Cj.  +  j  has  not  arrived  yet  when  leaves  the  sys¬ 
tem  (cf.  Figure  2.1b).  This  means  that  the  server  becomes  idle  at  the 
departure  of  Cj. ,  and  C^+  ^  arrives  at  the  system  after  a  delay  of  -  W^  - 
subsequent  to  the  departure  of  Cj,.  Let  Bj(t)  be  the  server  backlog  process 
of  class  I  M/G/l  system  when  Bj(0)  =  0.  When  C^+j  arrives  at  the  system, 
the  server  has  backlog  B^C-W^  -  While  the  server  works  on  his  back¬ 

log,  other  class  I  customers  may  arrive  who  also  precede  +  Hence,  when 
W’k  +  +  1  <  0,  the  waiting  time  of  C^  +  1  is  the  first  passage  time 

of  Bj(t)  from  its  initial  load  X  to  zero  where  X  =  Bj  (-Wj.  -  £^+j).  for 
notational  convenience,  we  define  Tqq  =  0. 

eff 

The  negative  support  of  the  variate  W^  +  £j.  +  i  =  W^  +  -  A^.  +  , 

is  contributed  only  by  which  is  exponentially  distributed.  Hence 

the  p.d.f.  of  W^  +  on  the  negative  real  axis  is  of  the  form 

*1  ix 

PkAIIe  U(‘x)  with  Pk  =  P^"k  +  ^k+1  "  where  Ufx-  =  1,  x  5  0,  and 

U(x)  =  0,  x  <  0.  Hence  when  W^  +  +  ^  <  0,  one  has 


-9- 


(2.7)  X  i  B j  ( -IXj.  -  htl)  i  BjCji-E)  ,  ht  ♦  ttl,  <  0  , 

where  E  is  the  exponential  variate  of  mean  one.  From  this  we  have  the 
following  theorem. 

Theorem  2. 1 

Let  be  the  waiting  time  before  first  entry  into  service  of  the 
k-th  class  II  customer.  Then 


(2.8)  »k„  =  < 


"k  *  5k.l  if  "k  *  Vl  *  °  ’ 


xo 


if  Kk  +  5k+l  <  0  ’ 


C 

where  +  j  =  S^.  -  A^+1  and  T^Q  is  the  first  passage  time  of  the  server 

backlog  process  B^(t)  for  class  I  M/G/l  system  from  X  to  rero.  The  distri¬ 
bution  of  X  is  given  by  (2.7). 

The  process  in  (2.8)  may  be  called  a  modified  Lindley  process  in 

that  the  homogeneous  process  +  1  is  modified  by  a  retaining  boundary 

at  the  origin  with  replacement  at  independently  chosen  T..^.  We  note  that 

-wX 

the  Laplace  transform  SY(w)  =  E[e  ]  of  X  in  (2.7)  can  be  written  as 

A 


-wBj(  E)  «  -X  t  -wB  (t) 

(2.9)  Bv(w)  =  E[e  ]  =  X  /  e  E[e  ]dt 

A  ii0 


“  _  t  -wBj  (t ) 

The  double  Laplace  transform  /  e  E[e  ]dt  is  given  [14,  22]  by 

0 


(2.10) 


■st  -wBI(t) 
StE[e  1  ]dt 


1  -  WE j (S) 

s  +  X'  {1  -  Oj  (w) }  -  w 


-10- 


Here  eT(s)  =  /  e"stE7(t)dt  where  E  (t)  =  P[B.(t)  =  0].  From  Equations  (2.91 
0  11 
and  (2.10),  one  then  concludes  that 


*TT{ 1  -  we  (A  ) } 

(2.11)  B..(w)  =  ^ - 

X  +  1  -  Clj  (w)  J  -  K 


The  following  theorem  can  now  be  readily  shown. 

Theorem  2 . 2 

-WTX0 

Let  oXQ(w)  =  E[e  ]  where  TX(7  is  given  in  Theorem  2.1. 

°X0M  “  '  w))[l  -  e^AjjHw  +  Aj  -  AjOgpj(w)}] 


Proof 


As  shown  in  (2.1)  of  [14],  one  has 


°X0(W)  6X(W  +  XI  '  XI°BPI^U^ 

An[l  -  eI(AII){w  +  A r  -  AjCgpjfw))] 

~  \u  ♦  Ajd  -  Oj  (w  ♦  Xj  -  AjOgpjfw))}  -  {*♦  Aj  -  Aj-bpjI'O 

The  theorem  then  follows  from  the  Takdcs  equation  (1.3).  0 

One  easily  finds  that  oXQ(w)  -*■  AjjEjfAjj)  as  w  -*  +»  So  that  Txo  has  mass 

^IIeI (^II)  at  the  ori8iru  Let  rX0(x)  be  probability  density  of  Txf)  on 

00 

the  positive  real  axis  and  define  yxq(w)  =  /  e_V,Xrxo (x)dx .  One  then  has 
YX0(”)  =  a  (k)  -  Ajjej(Ajj)  and  the  next  corollary  is  immediate  from 
Theorem  2.2. 

Corollary  2.3 

(a)  P[TX0  =  0]  =  P  [X  =  0]  =  AneI(AII)  . 


(fc)  YX0M  =  A 


n  -  w  f1  +  XIEIfXI 


l)cBpi(w)l 
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As  we  will  see  in  Section  4,  Corollary  2.3  plays  a  key  role  for  evaluat¬ 
ing  the  distribution  of  W'^  recursively  via  the  Laguerre  transform  method. 

The  numerical  value  of  is  therefore  needed. 

Theorem  2.4 

Let  h(s,w)  =  s  +  X  { 1  -  a^(w)}  -  w.  For  each  s  >  0,  h (s ,wl  is  strictly 
monotone  decreasing  in  w,  w  >  0,  and  has  a  unique  zero  at  w^fs)  =  l/ej(s)  >  d. 
Proof 

We  note  that  h(s,w)  =  (s  +  XT  -  w)  { 1  -  - r -  aT  (V)  / .  Then  from 

v  I  s  +  Aj  -  w  I 

Rouche's  theorem,  h(s,w)  has  only  one  zero.  This  unique  zero  is  attained 
on  the  positive  real  axis  when  s  >  0,  since 


^h(s,w)  =  -1 


*  -(1  - 


dt) 


<  0 


and  h(s,0)  =  s  and  h(s,w)  -+•  as  w  -*■  Hence  h(s,w)  attains  zero  at 
wQ(s)  >  0.  It  is  known  that  e^(s)  =  [s  +  A j  -  AjCgpjCs)]  .  One  then 
sees,  from  the  Takdcs  equation  (1.5),  that 


h(s, 


Gj  (S) 


5  *  h  • 


(5) 


1 

Cj(s) 


=  s  +  A. 


-  XIaBPI($)  -  S 


-  XI  + 


XI°BPI  ^  ~  ° 


Hence  w.(s)  =  — 
0  e  j 

Since  h(Ajj,  w) 
icallv  evaluate 


- — p,  proving  the  theorem.  C 
is  strictly  monotone  decreasing  in  w, 
Cj  (Aj  j)  =  1  /wQ (x j  r ")  straightforwardly. 


w 


^  0, 


one  can  numcr- 


Remark  2.5 


Ke  note  that  the  numerator  of  (2.10)  vanishes  at  w- (s)  =  — i-r-.  Hence 
“  _st  -KBj(t)  0  eI(sl 

I  e  E[e  ] dt  is  regular  for  Re(w)  >  0  when  s  >  0,  and  the  transfer;". 

0 

in  (2.10)  has  only  nonnegative  support.  This,  in  turn,  implies  that  ■,  (wl 

in  Corollary  2.3  has  nonnegative  support  only,  as  required. 

It  has  been  shown  (cf.  [10],  p.  235,  item  (d)}  that  the  ordinary 

Lindlev  waiting  time  process  is  stochastically  monotone.  In  particular, 

for  the  sequence  of  the  ordinarv  waiting  times  wK  one  has  h'J'  ,  h'!' . 

k  k+ 1  k 

k  =  0,1,2,...,  i.e.,  P[l\'k+j  >  x]  >  P[h'k  >  x]  ,  when  IVP  =  0.  The  waiting 
times  may  then  be  said  to  be  sequentially  monotone.  The  modified  Lindley 
process  (2.8)  is  also  sequentially  monotone  when  =  0,  as  we  prove  next. 
Theorem  2.6 


Let  WR  be  defined  as  in  (2.8)  with  =  0.  Then  are  sequentially 
monotone,  i.e.,  h'k+1>  h'k>  k  =  0,1,...  . 

Proof 


It  is  clear  that  IV  ^  >  IV  Q  =  0.  Suppose  >-  Kj.  ^ .  Let  =  [K^  +  £ 
where  [X]  =  max{0,X}  and  define  the  survival  functions  F..  ,  (x)  =  P[K,  >  x] 

I*  )  K  K 

and  fV  k^  =  P^Vk  >  xl*  "e  note  that  "k  ^  "k  1  imPlies  Vk  ^  Vk  1  ^cf' 
[10]).  From  (2.8)  one  sees  that 


(2.12)  F^k.jOO  -  Fv>i(x)  ‘  fVik(o.)-SX0(x)  .  x>0. 


Here  R^0(x)  =  P[T^q  >  x]  is  the  survival  function  of  the  replacement  distri 
bution  for  the  modified  Lindley  process.  From  the  induction  hypothesis, 
one  has  Fy  k(x)  ^  Fy  k  l(x),  x  >  0,  so  that  F^  kj.](x)  >  FK  k f x)  ,  x  >  0, 
proving  the  theorem.  □ 


Corol larv  2. 7 


(a)  Let  Ej.  =  PH'),  =  0] .  Then  =  F^,  j.  (0+)  •  Aj  jC  j  ( X  j  j  V  is  monotonical  ly 
decreasing  in  k. 

(b)  E[WjJ  is  monotonical ly  increasing  in  k. 

Proof 

Part  (a)  follows  from  Corollary  2.3,  Theorem  2.6,  and  (2.12).  Part 
(b)  is  immediate  from  Theorem  2.6.  C 
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customers;  its  heavy  traffic  approximation 

As  shown  in  Section  2,  the  sequence  of  time  spent  before  first 
entry  into  service  of  the  k-th  class  II  customer  is  a  modified  Lindley 
process  with  replacement.  For  the  ordinary  Lindley  process,  when  the 
virtual  value  1V^  +  is  negative,  replacement  is  at  zero.  For  our 

case  replacement  distribution  has  support  on  the  positive  continuum  as 
well  as  at  zero  and  the  standard  Kiener-Hopf  methods  must  be  altered. 

The  compensation  method  introduced  earlier  in  [6,  7]  and  presented  in 
simplified  form  in  [2],  provides  a  quick  analysis.  One  sees  that  the 
required  c.d.f.  F^.jj(x)  of  at  ergodicity  is  given  by 

00 

(3.1)  Fwn(x)  =  /  C(x-y)dc"(y)  . 

.00 

H 

Here  G^(x)  is  the  ergodic  green  distribution  of  the  underlying  homogeneou 
H  H 

process  K^  +  1  =  +  and  C(x)  is  the  c.d.f.  of  the  compensation.  In 

transform  notation  (3.1)  becomes 

(3.2)  4>WI  j  (w)  =  X(w)Y”(w) 


where 


(3.3) 


H,  •. 

Y„(”)  = 


II 


xii  •  “ 


eff .  . 

a  j  j  (w) 


and 


II 


X(„)  =  K[0SO(V)  - 


(3.4) 
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In  (3.5)  and  (3.4),  (w)  and  o^0(w)  are  described  in  (1.2)  and 

Theorem  2.2,  respectively.  K  is  a  normalization  constant. 

Consider  next  the  M/G/l  system  for  which  A  =  XT  +  A  and 
\  >  ’  bill 

AI  AI I 

Ac(x)  =  -r— A.(x)  +  t —  Att(x)  describing  the  service  stream  seen  bv 

b  As  I  As  11 

the  server  when  the  two  classes  are  given  equal  priority.  If  the 
service  discipline  is  FIFO,  the  Pollaczek-Khintchine  transform  for 
the  stationary  waiting  time  is 


1  -  P« 


(3.5)  WPKS^W^  * 


l-o 


1  -  as(w) 
S  wE  [T “] 


The  statement  in  (1.6)  may  now  be  given  formally. 

Theorem  3. 1 

Let  os  =  A<.E [Tg]  =  pj  +  <  1.  Then  the  ergodic  distribution  of 

the  time  until  first  entry  into  service  for  class  II  customers  with  FIFO 
discipline  has  the  transform 


(3.6)  ♦Kn(»0  s  WPKS (w  +  XI  ‘  XIaBPI(w)) 


Proof 

Let  Z(w)  =  w  ♦  Aj  -  AjOgpjfw).  From  Theorem  2.2  and  (5.4),  one  has 

X(w)  =  -KjZ(w)/ (A j j  -  w).  From  (1.2)  and  (3.3),  one  also  has 
X 

vV)  =  fl  -  t - — —  •  aTT(I(w))]"  .  Hence  from  (5.2),  one  sees  that 

»  v  1  A  -  -  -  W  II 


l_j  K  j  ^  ( u  1 

*IVII(")  =  =  w  -  Xn(l  -  a  j  T  ( Z  ( w ) ) } 
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Since  w  =  Z(w)  -  X^U  -  °gpj(K)}  =  "('0  -  Xj(l  -  oij(Z(w))},  one  then 
has 


KjZCw) 

♦WI(W)  =  Z(w)  -  XjO'-  aT  (Z(w) ) }  -  XTT{  1  -  aTT(Z(w))}  ‘ 


I 


II' 


II 


Dividing  both  the  numerator  and  the  denominator  by  Z(w),  the  above  equa¬ 
tion  leads  to 


Ww)  = 


1  -  os(Z(w)) 
1  "  PS  Z(w)EtTsi 


I  ^  1 1 

where  oig(wj  =  j-  ctj(w)  +  j —  djj(w).  From  6^II(0+)  =  1,  one  easily  finds 
S  S 

that  Kj  «  1  •  ps  and  the  theorem  follows  from  (3.5).  C 

From  Theorem  5.1  and  (2.5),  the  transform  of  the  time  in  system  of  class  II 
customers  at  ergodicity  is  immediate.  One  has: 

Theorem  3.2 

eff 

Let  Pg  <  1  and  let  =  fcjj  +  Sj ^  be  the  total  time  in  system  of 

jlw 

class  II  customers  at  ergodicity.  If  ^^(w)  =  E[e  ],  then 


(3.7)  j  j  ( w )  =  uip  (w  +  X  j  -  X  jO  gp  j(w))oijj(u  +  Xj  -  X  j  Ugp  j  (w ) ) 

The  transform  results  given  in  (1.2),  (1.3),  (3.5),  Theorem  5.1  and  Theorem 
3.2  lead  to  the  following  relations  between  the  first  moment  of  various 


variates  of  interest. 
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(3.8a)  E[TS]  =  EfTj]  +  E[Tn] 


E[T  ] 

<5-8b>  E<W  ■  1-^7  • 


fHnl 

1  '  Pt 


|3.8e)  EtK„]  =  E[KpKS](l  *  ^E[TBpI]>  = 


eff  EfKPKS]  +  E[TII1 
(3.8f)  E[Un]  =  E[Wn]  ♦  E (Tj j  ]  =  - - - ~ 


The  heavy  traffic  limit  theorem  for  class  II  customers  can  now  be  given. 
Theorem  3.3 


Let  the  competing  service  times  Tj,  Tjj  have  finite  second  moments. 

Let  (Aj.,  Ajj.)  be  a  sequence  of  arrival  rates  for  which  Aj  =  K  A ^ j  , 

J  J  i  i 

K  >  0,  and  let  p<.  -*•  1-.  Then  both  K^/E[K^]  and  U^/EfU^]  converge  in 

distribution  to  the  exponential  variate  with  mean  one. 


Proof 

When  Aj  and  A^  are  in  fixed  ratio,  the  system  service  time  T§  has 

j  j 

the  transform  oig(w)  =  (A^  /A^  )a^ (w)  ♦  (A^  /A^  )a^(w)  so  that  the  distri- 

i  j  i  i 

bution  of  Tg  does  not  change  with  j.  From  Theorem  3.1,  (5.5)  and  (3.Se), 


(3.8c) 


EUiild  .  x,E[TBpI]) 


(3. 8d) 


E[KPkS] 


°S  E^ 

2(1  -  P  )  '  E[T  ]  ‘ 


one  has 


-IS- 


(5-9>  wf->  ■ - 7 — \ — : - 

11  ‘■»s“s'f'l*rrrVilri1> 

II  1  II 


where  <j>^(w)  =  {1  -  <t>^(w)  }/wE[X]  and  E[Wjj]  =  Wjj-  Further,  from  (5.8) 


and  (3.9), 


(3.!°)  *KII<r->  ■  (>  -  os”'  -  -  V1  -  Vi(r-’liir'  - 

'll  PKS  'll 

* 

One  may  then  employ  the  Taylor  expansion  with  remainder  of  both  a  (w)  and 

★ 

Ogpifw)  out  to  the  linear  term  in  w,  and  proceed  classically  invoking  the 
continuity  theorem  for  characteristic  functions  to  find 

'3-u>  Wr-’-Tk'  as  °s  *  '•  • 

"il 

demonstrating  the  convergence  of  Kjj/E[Kjj],  The  convergence  of  Ujj/EfUjj] 
is  immediate  since  the  effective  service  time  is  bounded  stochastically 


from  above.  □ 


Remark  3.4 


One  could  also  inquire  about  the  limiting  behavior  when  Aj  and  A^ 

i  i 

are  not  in  fixed  ratio.  For  application  of  heavy  traffic  approximation, 
however,  one  has  specified  values  of  A^  and  A^  with  os  =  1-e  for  e  >  0, 
small.  The  approach  path  of  A^  and  A^  to  the  values  A^  and  A^  is  then 

1  i  1  i 

irrelevant. 
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§4.  Evaluation  of  the  distributions  of  interest  via  the  Laguerre  transform 
In  this  section  the  distributions  of  interest  described  in  the  pre¬ 
vious  sections  are  evaluated  numerically  via  the  Laguerre  transform.  The 
reader  is  referred  to  [11,12,21]  for  the  underlying  theory  of  the  Laguerre 
transform.  All  figures  are  given  at  the  end  of  this  section.  The  follow¬ 
ing  example  is  considered: 

System  I:  High  priority  class 

*j  =  \  ,  Aj(x)  =  P[Tj  >  x]  =  |[e'"2X  +  e'x  +  e-;,x] 

E[Tj]  »  ~  ,  Pj  -  AjEfTj]  =  0.278  . 

System  II:  Low  priority  class 

An  *  j  ,  Ajj(x)  =  P[Tn  >  x]  =  j[e"X  +  e'“X] 

E  [Ti  i  ]  =  |  ,  Pn  =  AjjEfTjj]  =  0.375  . 

Total  system 

XS  =  Xl  +  AII  =  4  ’  ^S(X^  =  PtTS  >  =  J  *1^  +  3  ^11^ 

E[TS]  -  i  E[Tj]  ♦  |  E[Tn]  =  0.870  ,  Pg  =  Oj  ♦  Pn  «  0.653  . 

(A]  Ergodic  distributions  of  and 

From  the  transform  results  (1.2),  (1.3)  and  (5.5),  one  obtains 
easily  the  corresponding  (generalized)  probability  densities  as  given 


below. 
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(4.1) 


00  -X  x  (X  x)n 

s r>n t  (x)  =  l  e  -n7T)!  ain*1  (x) 


BP  I 


n=0 


(4.2) 


££  00  -X  X  (X  \) 

eff,  r  ,  I  I 

stt  M  =  /.  (e 


(n) 


n=0 


n!  '  an(x)}*sBpi(x) 


(4.3) 


fPKS^ 


=  (1  - 


cs)6(x)  +  fpKS(x) 


where 


(4.4) 


fPKSW 


=  (1  - 


PS'1  DS  {E(TS] 


Here  a^(x)  is  the  n-fold  convolution  of  a(x)  with  itself.  The  asterisk 
also  denotes  convolution  and  6(x)  is  the  delta  function. 

Similarly  from  Theorem  3.1  and  Theorem  3.2,  the  probability  densities 
fWn(x)  and  fyjj(x)  take  the  forms 

(4.5)  fKn(x)  =  (1  -  ps)6(x)  +  f*n(x) 


where 


(4.6) 

and 

(4.7) 


fKII W 


l  te 

n=0 


-XjX  (Xjx) 


n  i 


^PKS  ^  ^*SBPI  ^ 


fUII(x) 


,.  ..  eff,  .  e*  r  \  +  eff, 

=  (1  -  os)sn  (x)  +  f IV! I  x^ * s 1 1  (x^ 


The  Laguerre  transform  enables  one  to  evaluate  systematically  all  of 
these  probability  densities  which  heretofore  have  been  behind  "the 
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Laplacian  curtain".  The  Fourier-Laguerre  sharp  coefficients  of  a^fx) 
and  a^(x)  are  readily  obtained  analytically.  Using  the  relevant  opera¬ 
tional  properties  of  the  transform,  Equations  (4.1)  through  (4.7)  lead 
to  the  coefficients  of  each  density.  They,  in  turn,  can  be  converted  to 
the  coefficients  of  the  corresponding  survival  functions,  thereby  bypass¬ 
ing  numerical  integration.  The  final  inversion  from  the  Laguerre  coeffi¬ 
cients  to  function  values  is  straightforward.  In  Figure  4.1,  the  survival 
functions  Sgpi(x)  =  P[Tgpi  >  x]  and  S^f(x)  =  P f S ^ ^  >  x]  are  plotted. 
Figure  4.2  depicts  the  survival  functions  =  P[Kp^<-  >  x]  , 

FWII(x)  =  P^VII  >  and  PUII^  =  P^UII  >  x^'  lVe  note  that  both  kpks 

and  Kjj  have  mass  1  -  p<.  at  the  origin. 

(B)  Modified  Lindley  process  with  replacement 

It  has  been  seen  in  Theorem  2.1  that  the  waiting  time  before  first 
entry  into  service,  W^,  of  the  k-th  class  II  customer  follows  the  modified 
Lindley  process  given  in  (2.8),  i.e.. 


(4.8)  Kk+1 


"k  +  Ck+1  lf  Wk  +  5k+l  "  0  ‘ 


if  Wk  +  5k+l  <  0  ’ 


-Tv 

V, 


where  -  A^+j  and  the  transform  o^  (w)  =  E[e  0  ]  is  given 

in  Theorem  2.2.  As  shown  in  (2.3),  the  total  time  spend  in  the  system 
by  the  k-th  class  II  customer  is  then  given  by 

(4.9)  U„  -  Wk  .  s'ff  . 
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We  now  show  how  these  transient  distributions  can  be  evaluated  via  the 
Laguerre  transform. 

Let  a(x)  be  the  p.d.f.  of  the  i.i.d.  random  variates  £ The  Laguerre 

t  CO 

sharp  coefficients  (a^)  x  of  a(x)  are  easily  obtained  from  those  corre- 

£  f  _ 

sponding  to  Sr  and  A  j.  From  Corollary  2.3,  the  variate  has  mass 

R0  =  at  the  origin,  which  can  be  calculated  using  Theorem  2.4. 

In  our  example,  RQ  =  0.848.  The  probability  density  r^  (x)  of  T  on  the 

positive  real  axis  has  the  Laplace  transform  y  (w)  given  in  Corollarv  2.3. 

r  CO 

The  Laguerre  sharp  coefficients  Crn)Q  of  r^Q(x)  are  tEen  obtained,  using 

Corollary  2.3(b),  from  those  corresponding  to  T^  and  Ar.  bet 

Er  =  P[Wr  =  0]  and  let  fk(x)  be  the  probability  density  of  Kr  on  (O,^)  so 
00 

that  E,  +  /  f,(x)dx  =  1.  Assuming  that  E,  and  the  Laguerre  sharp  coeffi- 

#0  oo 

cients  (f'(k))_  of  f.  (x)  are  known,  we  next  establish  an  algorithm  for 
n  0  k 

obtaining  Er+1  and  (f*(k+l))“  in  terms  of  (a')*^,  RQ,  Ek  and 

#  ao 

(fnOO)0. 

U 

Let  f"+1  (x)  be  the  probability  density  of  Wr  +  ^k+1'  tEen  Eia? 

(4.10)  fR+J(x)  =  EkaW  4  fk(x)*a(x)  ,  -»  <  x  <  »  . 

QO  J) 

The  associated  Laguerre  sharp  coefficients  (t  (k+1))  of  f!  . (x)  are  then 
p  '  n  k+1 

given  by 

(4.11)  f  (k  +  1)  =  E,  a  +  7  a  f  (k)  ,  -®  <  n  <  °° 

’  n  v  ’  k  n  n-m  mk 

m=0 


Let  fk+1(x)  =  fk+j(x)U(x)  where  U(x)  =  0,  x  <  0  and  U(x)  =  1,  x 

+  rt  oq  ^ 

coefficients  (f  (k+l))„  of  f,  , (x)  are  found  from 
n  0  k+1 


>  0.  The 
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(4.12) 


f+0Vn  =  -  I  fJJ^Ck+l)  ;  f**  (k+l)  =  fH=(k+l)  ,  n  >  1 

n=l  ''  11 


Let  pk.l  *  pl*l  *  «k.l  •=  °I-  The"  Pk*l  ■  I  fk.l<x)dx  -  1  -  / 


so  that 


(4-i3)  p"  ,  =  i  .  :  l  . 


k+l 


n=0 


2n+ 1 


From  (4.8),  one  finally  has 


(4.14) 


:k+l  ~  Pk+1R0 


and 


(4.15)  f*(k+l)  =  fn"(k+1)  +  P^g*  ,  n  >  0  . 

Equations  (4.11)  through  (4.15)  enable  one  to  calculate  E,  ,  and  (f^ (k+ 

k+l  n 

recursively  for  k  =  0,1,2,...,  starting  with  Kp  =  0  (i.e.,  Ep  =  1,  fy.u)  = 

#  *  oo 

and  fn(0)  =  0,  n  >  1).  The  coefficients  Cfyn (k ) ) 0  corresponding  to  are 
given  from  (4.10)  by 


t.  ,,,  c  eff*  r  *  ,  eff* 

(4.16)  f  (k)  =  E.  s  +  >  f  (k)s 

Un  k  n  Ln  n-nr  ’  m 

m=0 

where  (s®^  )p  are  the  Laguerre  sharp  coefficients  of 

In  Figure  4.3,  the  survival  functions  F^(x)  =  P[l\^  >  x]  are  plotted 
for  k  =  1,2,3,4,5,10,20,30,40,50  and  0  <  x  <  10.  The  absolute  difference 


between  F(^.,.p(x)  and  its  crgodic  survival  function  FIt(x 


I  I 


wn 


fvid- 


,-6 


obtained  in  (A)  is  bounded  by  10*  for  0  <  x  <  10,  usinc  the  first  130 
Laguerre  coefficients.  This  assures  numerical  stability  and  accuracy 
of  the  Laguerre  transform  procedure.  Ke  note  that  the  stochastic  mono- 
tonicity  of  in  k  given  in  Theorem  2.6  can  be  observed  in  Figure  J.3. 
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Figure  4.4  and  Figure  4.5  show  the  convergence  of  Ej_  to  1  -  o<~  and  that 
of  E[Wj.]  to  E[Wjj]  as  k  -*■  ®,  respectively.  Both  E^  and'E[h'^]  are  calcu¬ 
lated  using  the  Laguerre  sharp  coefficients.  It  has  been  shown  [15,  21] 
that  the  Laguerre  sharp  norm  defined  by 


(4.17) 


provides  a  distance  between  any  two  distributions.  In  Figure  4.6  this 
Laguerre  sharp  norm  distance  between  and  for  1  <  k  £  50  is  exhibited 
thereby  quantifying  the  rate  of  approach  to  ergodicity.  These  distances 
also  provide  convenient  stopping  criterion  for  the  computation.  One  can 
see  that  for  k  >  25  the  distance  is  bounded  by  0.01.  Finally,  in  Figure 


4.7  the  survival  functions  =  P [U^  >  x]  are  plotted  for  k  =  1,2, 5, 4, 

5,10,20,30,40,50  and  0  £  x  £  10.  All  computations  were  done  on  a  PEC  10 


computer  in  a  timesharing  mode  using  APL  as  the  programming  language. 


Relevant  formulae  were  usually  coded  in  a  straightforward  way  using  the 


first  150  Laguerre  coefficients.  The  tesults  displayed  here  were  typically 


obtained  with  CPU  time  in  seconds. 
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